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Abstract. We calculate a mass of a quintessence particle of order 10 -5 eV and we 
find several solutions for quintessence field equation. We consider also a quintessence 
^ | speed of sound in several schemes and quintessence fluctuations. 
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In this paper we consider several consequences of the Nonsymmetric Kaluza-Klein 
' (Jordan-Thiry) Theory. We consider a value of a mass of quintessence particle, sev- 

Q ■ eral interesting relations among energy scales, radiation density in the second de Sitter 

phase. We find a spatial dependence of a quintessence field (and an effective grav- 
£f) ' itational constant C7 e fr) in a case of spherical-statical symmetry, cylindrical-statical 

symmetry, flat-static symmetry. We find a time dependences of a quintessence field 
(with no spatial dependence). We get a solution for a quintessence field (a travelling 
wave) and two-dimensional wave solution applying those solutions for G e s . We propose 
some kind of statistical approach to our results. We calculate a speed of sound in a 
quintessence. We consider also fluctuations of a quintessence caused by gravitational 
waves perturbations. 

Let us consider a value of mass of a quintessence particle (a scalar particle) 
(see Ref. [1]) obtained from Nonsymmetric Kaluza-Klein (Jordan-Thiry) Theory (see 
^ ' Ref. [2]): 

^( n \ n/ \ , 1/2 f\^ n/i 
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(see Eq. (79) from Ref. [1]). 

The value of this mass has been obtained by this particle during the second de 
Sitter phase. Moreover during our contemporary epoch it is the same. The parameters 
^ n, 7 and are defined in Refs [1], [2]. During the second de Sitter phase a cosmological 

constant has been calculated in Ref. [2] and one finds 

2 2|7|( n+2 )/ 2 n™/ 2 
A c0 = 6i/ 1 = pn/2{n + 2){n+2)/2 . (2) 

The value of a cosmological constant remains the same during our contemporary epoch. 
Hi is a Hubble constant during the second de Sitter phase (see Eq. (8) of Ref. [1]). 
Thus for our contemporary epoch one gets 

m o = 2 V n ( n + 2 ) \AcO • (3) 
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In this way we can connect the value of a cosmological constant of our contemporary 
epoch to the value of a mass of a quintessence particle. From recent observational data 
we get 

1 



A c0 = /l = 10- 52 ^. (4) 



ur 

Moreover in order to get a correct dimension for a mass we should add a factor with h 
and c (now we abandon the system of units with h = c = 1). 
One gets 

1 H /— r 

m o = „- yn(n + 2) V A c0 . (3a) 
^ c 

Using the value of A c q one finally gets 



m 



~ V«(« + 2) -0.17- 10" 39 g (5) 



or 



m ~ + 2) • 0.95 • 10 -5 eV. (5a) 

For example, if we take n = 14(= dimG2), one finally gets 

m ~ 14.2 • 10" 5 eV. (6) 

This value is bigger than that considered by different authors. Moreover, still 
sufficiently small. The particle interacts only gravitationally and because of this it is 
undetectable by using known experimental methods. Taking a density of dark energy 
as 0.7 of a critical density, 

Pc = 1.88ft 2 • 10" 29 -^, (7) 

cm. 6 

one gets a number of quintessence particles per unit volume 

n= - k2 .l.31-10 10 -L ( 8 ) 

y/n(n + 2) cm 3 

where h is a dimensionless Hubble constant 0.7 < h < 1. Taking n = 14 and h = 0.7 
one finally gets 

n = 4 • 10 8 -Xt (8a) 
cm 3 

which is many orders of magnitude smaller than Loschmidt number. Thus a gas of 
quintessence particles is not so dense from the point of view of our earth conditions. 
However, if this number of particles per unit volume is considered in a container of size 
200 Mpc, the gas can be considered as extremely dense. 
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In order to settle — is this gas dense or not — we should calculate a mean scattering 
length. The scattering cross-section for a quintessence particle 



A mean scattering length 



a = — = 10 52 m 2 . 



I = — 

an 



where n is a number of quintessence particles per unit volume (Eq. (8a)). 
One gets 

/ = l(n 60 m. 



(*) 



(**) 



(***) 



It means that a gas of quintessence particles is extremely dense (if we apply the Knudsen 
criterion — a gas is dense if I <C L, where L is the size of the container) even in the 
Solar System. 

Let us consider the Eq. (Ill) from Ref. [3] which connects several scales of energy 
and gives an account of the smallness of gravitational interactions in our theory. We 
rewrite this equation in the form 



( m P' ^ ( m A 



(ra+2)(m+2)/2 



m E w/ \(n + 2)/3 / 

where m^w is an electro-weak interactions energy scale. Taking 

= m EW 
m = 2 (M = s 2 ) 
n=U(H = G2) 



(9) 



one gets 



where 



21 



and eventually one gets 



K 



( m p' ^ 



1/24 



Taking ibew — 80GeV and m p \ ~ 2.4 • 10 18 GeV one gets 

k ~ 6.19 



(10) 
(11) 
(12) 

(13) 



3 



which is very reasonable for it establishes a relation between two cosmological terms 
7 and j3 as of the same order. Simultaneously this is a consistency condition for our 
model with energy scales (the 6-dimensional Planck's mass is equal to uiew)- In this 
way we can calculate a mass of a quintessence particle for our contemporary epoch 
from Eq. (1) 

m = 2- 10 15 eV- |£| 1/2 . (14) 
This gives us an estimation for \P] and R{r): 

R(T) = -( -^W-V |P| (15) 
k \m pl a ein J 

where a em = is a fine coupling constant. Using Eq. (13) and values of m p \ and 
ttt-ew one gets 

R(r) = 24.75- 10- 31 |P|. (16) 
From Eq. (14) and Eq. (6) one gets 

\P\ = 5 ■ 10 -39 (17) 

and 

R(r) 9* 1.2 • 10" 68 . (18) 
Moreover in our simplified theory we have 

~ ~ 2(2/i 3 + 7 M 2 + 5^ + 20) 

R{F) = WTW ( } 

and [i should be very close to the root of the polynomial 

W{ii) = 2/x 3 + 111 2 + 5fi + 20. (20) 

From (18) and (19) one gets 

2^ + 7^ + 5^ + 20^ 1.3- 10" 66 . (21) 

Thus /i is very close to the 70-digit approximation of the root of the polynomial (20) 
(W(jl) = 0.1 • 10 -67 ). Due to this we can control the cosmological terms. 
In the case of \P\ we have the formula (28) from Ref. [1] and one gets 

~ ~ dP 

P((o + e) = P((o) + ^(to)e (22) 

£(Co = ±1-38 . . . ) = (23) 
dP 

^=(| Co | = 1.38...) = 25. (24) 
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Thus one gets from (17) and (23-24) 



e ~ 2 • 1(T 4U . (25) 

Thus we need an approximation of Co U P 40-digit arithmetics. We see that cosmological 
terms coming from the Nonsymmetric Kaluza-Klein (Jordan-Thiry) Theory are very 
small, but not zero, and that they are easily controllable by fi and ( parameters. 

Let us consider a self-interaction potential for a quintessence field for our contem- 
porary epoch (which is the same as for the second de Sitter phase). One gets from 
cosmological terms A c o(^) 

M*0 = ~ - |7l) e"* (26) 

rrf \ \l\ ( Al\ V I nm pl 

U{qo) = ~ TTi T~o\ ( exp / go 

2(n + 2) \(n + 2)0j \ 2 . 2it\~M\ 



x I n exp I ^— q I - (n + 2) I (27) 



or 



U(q ) = -tAcO exp -^^q nexp -fi^q, I - (n + 2) I (27a) 
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2^2n\M\ J \ ' \^2ir\M\ 



where A c o is a cosmological constant for our contemporary epoch and 

V = % + — mp ' go = #b + Pqo = % + <P- (28) 



2y2vr|M| 

For A c o is very small (lO -52 ^-), this interaction is very small. For n = 14(= dimH = 
dimG2) one gets 

U(qo) = -^A c0 exp f -^U/o j | 7«>xp | -pH=<,» | - 8 | . (29) 
V^2tt|M| 

Even in the potential U(qo) we have exponential terms, the strength of the inter- 
actions is negligible for small value of go- The interesting point in our theory is the 
effective gravitational constant (depending on a scalar field Let us describe it by a 
quintessence field go- One gets 

G^ = G e^ n+ ^. (30) 
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After some calculations one finds 

f (n + 2)P \ n+2 ( n + 2 Try , 
G e ff = G n — exp p g0 (31 

or 

G 4(n + 2) 2 /3 2 / n + 2 m pl 

Geff = T2 -2 ex P ~ H= 9 ° | • (31a) 

A c0 n y 1 \J2tt\M\ 

It is easy to see that the role of a Newton's constant Gn is played by 



G /2(n + 2)/3 



2 



a » = t\^r^ ] ■ (32) 



This is the reason that we put a constant Go in the formula (30). 
One gets 



G -G N X 2 c0 [ 2{n l 2)(3 ) • (33) 



For 



one gets 



or 



/3 = a> 2 1J R(r) = a 2 ^ (34) 



Gk _ G^/2&i±aa©y (35) 

(-Wpi) 2 In/ 



Go = G N [ |-^| • (36) 

l2(n + 2)i?(r) 1 ' 



Let us connect to Go a new Planck's mass m%. In terms of this mass and ordinary 
m p i mass one gets 

2(n + 2)fl(f ) a 2 
m P' " mpI n A^J (37) 

or 

m p i = m", ^3 ^- . (38) 

P 2(n + 2)R(r) al 

If we take our simplified model with R(r) given by Eq. (19) 

(X s ftein T37 ' ^ 
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we get 

o 32(2^ + 7^ + 5^ + 20) a em 

mpl = mpl WTW (39) 



or 



7(^ 2 + 4) 2 \V A ^pi 



G °- G H^wr^jJ l^rj • (40) 

Moreover we can connect unknown constant Go {m%) with an energy scale of electro- 
weak interactions tbew in a way suggested in Ref. [3] and developed here in this paper. 
It is enough to use the formula (10) and we get 

f7n\ 24 32(2// + 7/j 2 + 5/J + 20) a em 

mpl = mEW UJ rWTW W (41) 

If we take for A c o the known value of the cosmological constant (Eq. (4)) for l p \ = 
4 • 10~ 32 cm and k = 6.19 we get 



m° -8 19 10 132 V + V + 5/^ + 20 

m pl -8.12-10 m EW - ( 1-2) 



Moreover we see that /j must be very close to the root of the polynomial 2/j 3 + 7/z + 
5/i + 20 (see Ref. [3] ) and we get 

m°, = 2 ■ 10 130 (2^ 3 + 7/j 2 + 5/i + 20)m EW • (43) 

If we take 70-digit approximation of the root of the polynomial (20) and considered 
value of W(n) (Eq. (18)) obtained here, we get 

m° pl = 4.6 • 10 64 m EW (44) 

or 

m° = 4.6 • 10 64 m p i = 1.5 • 10 48 m pl = 3.2 • 10 43 g = 1.6 • 1O 1O M (45) 

p m p i 

where M Q = 1.99 • 10 33 g is a Solar mass. (This is a mass of our Galaxy.) 
Moreover the present critical density of matter in the Universe is 



ft ,„ = 2.775ft- X 10" IF *L_ (46) 



So we can find a volume containing m pl . One gets 

Pc,0 



y = — = 5- 10" 3 /i 2 (Mpc) 3 . (47) 



Thus the size of this volume is of order 



L ~ v^O.lTMpc. (48) 

Taking under consideration the fact that our calculations are quite rough (for 
example k could be bigger a little, or Wijl) a little bigger), we come to the conclusion 
that is of order of the mass of our visible Universe (a volume of size of 10 4 Mpc). 
In this way it seems natural to suppose that 

m° pl = M v (49) 

where My is the total mass of the visible Universe. However, this intriguing conjecture 
demands many investigations. 

Moreover it is in a real spirit of Dirac's large number hypothesis and can give a 
link between cosmology and fundamental interactions theory. Thus if we suppose (49) 
we get 

Go = G u = J L. (50) 
Let us come back to the equation for an effective gravitational constant (Eq. (32)): 

G e ff = G N exp I - n ^~ , mpl qo I • (51) 



qo — a quintessence scalar field — possesses a mass and because of this it has a finite 
range with Youkave type behaviour 

« ( R \ 
q = n exp (52) 



R L \ r 
where 

2 

r = — _ and a is a positive constant (53) 

v / n(n + 2)vA^ 

or taking for n = 14 and for A c o Eq. (4), 

r = 3- 10 23 m= 10 Mpc. (54) 

The formula (52) is valid for R ~ r$. Thus we cannot observe scalar gravitational 
radiation from closed binary sources and the quadrupole radiation formula is satisfied 
for a gravitational field. 

Moreover for R < r we should take under consideration a full self-interaction 
potential of a field qo (Eq. (29)). Moreover, because of the constant A c o in front of the 
formula (29) this is negligible. Because of this we can repeat some considerations from 
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the first point of Ref. [2] coming back to the field & and consider different sources of 
mass for this field due to interaction with the matter. In this way in both cases 



GW = G N . (55) 

However, we can expect some small effects on short distances. 

Let us notice that in the previous approximation we consider a weak field, it means 
| go I *C 1- This is different from small field. A small field considered below is such that 
qo is small in a usual sense. It is negative. Moreover it can be big in the sense of the 
absolute value. 

Let us calculate a radiation density in the moment when the second de Sitter phase 
starts. It means the radiation released after the phase transition. One gets 

Pr = Aci(#i) - A c0 (^ ) = A c0 (^o) (§| - l) ■ (56) 
Using Eqs (la) and (8a) from Ref. [2] one gets 



A c0 (^o 



n 2 



\P\ 2 + \P\\Jn 2 \P\ 2 + 4(n 2 - A)AR{r) + AAR{r){n + 2) 

2(n+2)/2|J>|(n+2)/2 n n/2 

(n-2)/2 



(57) 



x n 



(n\P\ + \Jn 2 \P\ 2 + 4(n 2 - 4)AR(r)^ 

/m~\ n m 2 ~ I p|(ra+2)/2 n/2 

M*o) = «n? = 2(^) ^ f—r 2 . 

V"W a s " (n + 2) (n + 2)/2 (R(r)\ 



(58) 



The density p r is of course an effective radiation density, because we adsorbe into p r a 
factor with an effective gravitational constant. 
In our simplified model we get 



A 



' 98g 2 (C, p) + (p 2 + 4) yjog ((, p) + 384fr(C P) + 32ff (C, p) 
g8((,p)(p 2 + 4) 2 -2^-7™ 

x (7g((, p) + v /49 ff 2 (C, / u) + 384/ l (C, / u)) ' (59) 
x (ln(|C| + v / C^ + T)+2C 2 + l) 6 -l^, 
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m pl J af 2 24 
v g 8 (C,/i)(M 2 + 4) 8 ( 6 °) 

X g . 

(2^3 + 7)U 2 + 5fi + 20)7 (in (|C| + v'C r +T) + 2C 2 + l) 

In the formulas (59-60) we should put = 7Bew and o? s = a em = ^?f- Let us notice 
that A is our contemporary cosmological constant given by (4). It gives a scale for p r . 

The functions g{C,n), h((,n), H(C,fi) are given by Eqs (40), (41a), (41b), (42) 
from Ref. [1]. The numerical factor in front of A can be calculated and we get 

A = 1.04- 10- 178 m| w a(C,/x) (61) 

where 

««, ,) = g8<C + 4>S , ■ (62) 

(2^3 + 7^2 + 5^ + 20)7 (in (|C| + a/C^TI) + 2C 2 + l) 

The calculated radiation density evolves in time according to the theory developed in 
Ref. [4]. 

In order to find some influence of qo (quintessence) field on the value of the effective 
gravitational constant we consider a field equation for the scalar q field in empty space. 
One gets 

d 2 qo 



. )r2 - V 2 (7o) + eaexp(n(3q ) (exp(2/fy ) - l) = (63) 



where 



_ A c0 n(n + 2)m 2 1 

a = = — — (64) 

16-7T • M 

P = ^= (65) 
2^2tt\M\ 

e = sgnM, e 2 = 1. (66) 

Let us consider a static, spherically symmetric case. In the spherical coordinates 
one gets 

° = hlk if IF ) " ^ ex P( n ^o) (exp(2^ ) - 1) (67) 

where qo = qo( r ) is a function of r only. In order to treat this equation it is easier to 
come back to the old variable ip = (3qo (see Eq. (28)). One gets 



^('• 2 ^) 7-"'' () - (cs) 
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We change the independent variable r into r 

^ (^) " eexp(n^) (exp(2^) - 1) = (69) 

where 

r = —2= t (70) 
V A c o 

and 

T n(n + 2)A c0 2 , , 

A c0 = = m„, . (71) 

8vrM pl 1 ; 

We consider Eq. (69) in two regions: 

1) for small fields ip, 

2) for large fields <p. 

In the first region we get 

1 d f T 2^\ + ~ e n, =Q (?2) 



r 2 dr V dr 
In the second region we get 

1 df 2 d^\_ Ie(n+2)v = Q 
t z dr \ dr J 

Let us notice that both equations have similar nature and can be reduced to the 
equation 

where in the first region e = 1, 

y = (75) 

x = y/nr, (76) 

and in the second region e = — 1, 

y = (n + 2V, (77) 
x = + 2 r, (78) 

We can transform (74) into 

x fy +2 ^- + ee xe y = (79) 
dx A dx 
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which is the celebrated Emden-Fowler equation known in the theory of gaseous spheres 
(see [5]). Let us notice that the first region (small fields) means large distances and 
the second region (large fields) means small distances. 

In this way we should consider Eq. (79) in the region of small and large x. In the 
case of ££ = 1 the equation (74) has an exact solution 



y = ln [-* 



j- 



(80) 



Let us apply this to both regions (remembering that e in both cases has a different 
sign). 

One gets in the first region 



<7o 



m p \n 



I 



ln 



2V2 



and 



G e ff = G 



N 



\ 



\ V n\ c0 
(n+2)/n 



2V2 



In the second region 



Qo = - 



2^2n\M\ 
m p i(n + 2) 



ln 



2V2 



and 



V(n+2)A c0 / 



G e ff = G 



N 



2V2 



V(n+2)A c0 / 



(81) 



(82) 



(83) 



(84) 



In this way we get an interesting prediction for the behaviour of the strength of gravi- 
tational interactions. In this very special solution G e g is going to zero if r — > and to 
infinity if r — ► 00. 

Let us come to Eq. (79) supposing ee = 1. Thus we get 



d2 y , n d y , v n 

x —, j + 2 — + xe y = 0. 

ax z ax 



(85) 



Using an exact solution (80) we write 



y = yi + : 



(86) 
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and consider Eq. (85) for large x. 

In this way we get an approximate solution (given by Chandrasekhar [5]) 

y = hi(Jp) + Acos^ln^ -2\nS, \A\ « 1, (87) 

where rj = ■=-, A and 5 are integration constants, 5 > 0. In this way we get in the first 
region 

GefT = G eS ■ exp ^- cos In 77^ ^ (88a) 
where yl is a constant (|A| <C 1) and 



V = ^vr r > (88b) 



G e ff is given by the formula (82). 
In the second region 



GW = G eff • exp ^- cos In 77^ ^ (88c) 

where A is a constant (|A| <C 1), G e ff is given by the formula (84) and 

^(n + 2)X c0 5 



2^/2 



(88d) 



In this way we have very interesting non-Newtonian behaviour of G e ff for large dis- 
tances. Let us notice that the length scale is completely arbitrary, because it is given 
by an integration constant 5. 

Let us consider Eq. (63) in Cartesian coordinates supposing flat symmetry for a 
quintessence field qo = qo(z,t) (nonstatic). One gets 

~W ~ ~dz^J ~ ? « ex P( n ^°) ( ex P( 2 ^o) - 1) = 0. (89) 
Let us change dependent and independent variables to £, 77, ip: 

z = ^=i (90) 
V X c o 

t = ^=r) (91) 

<P = 0qo- (92) 
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One gets 



(93) 



Eq. (93) is an equation for flat scalar (quintessence) waves in our theory. Let us consider 
it for large and small field (p (as before). 
In this way one gets the equation 



d 2 y d 2 y 
dT2 ~ dx 2 



+ eee^ 



0, 



where in the region of small field ip we have e = 1 and 



(94) 



y = nip 
X = y/n£ 
T = y/nrj 



(95) 
(96) 
(97) 



and in the region of large field (p, e = — 1 and 



y = (n + 2)ip 
x = Vn + 2^ 
T = Vn + 2r). 



(98) 
(99) 
(100) 



Eq. (94) is the famous Liouville equation which can be transformed via a Backlund 
transformation into a two-dimensional wave equation and afterwards solved exactly. 
The general solution depends on two arbitrary functions / and g of one variable, 
sufficiently regular. It is possible to consider several problems for this equation: Cauchy 
initial problem, Darboux problem and Goursat problem. 

The general solution of (94) looks like (ee = —1) 



y(T,x) =ln 



2g'(x-T)f'(x + T) 
(g( x -T) + f(x + T)y 



where g' and /' are derivatives of g and /. 

Thus one gets in the first region (small field) 



(101) 



2v^M 

q {t,z) = In 

mpin 



2^2 



(102) 
(103) 
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and 

G e ff = Gn I 

In the second region (large field) 



(n+2)/n 



, , 2V2vrM 
q (t, z) = ; — • In 



(104) 



m p i(n + 2) 



6 = 



2^ 



and 



G e ff = Gn ■ 



^(n + 2)A c0 



(105) 
(106) 

(107) 



In this way we get a spatio-temporal pattern of changing the effective gravitational 
constant for small and large field regions. In both cases we have ee = —1. However, in 
the small field region we have e = 1 and because of this e = — 1 . In the case of large 
field region e = — 1 and e = 1. In order to be in line with our assumptions we should 
consider in the first case such functions / and g that the expression in (105) is small 
and for the second case vice versa. 

Let us consider Eq. (63) in cylindrical coordinates supposing cylindrical symmetry 
for the field (?o> <Zo = Qo(p)- One gets 



dqp 
dp 



— ea exp(nf3q ) (exp(2/3g ) — l) = 0. 



This equation can be transformed into 



1 _d_ 

t dr 



r^\- eexp(n^) (e 2 ^ - l) = 0, 



A. 



c0 



As usual we consider Eq. (109) in two regions for small and large fields. 
In the first region 



r dr 



dip 



+ ee n <? = 0. 



In the second region we get 



1 d 

t dr 



dip 



ee 



(108) 

(109) 
(110) 

(111) 
(112) 
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Both equations can be reduced to the equation 



x dx \ dx 
where in the first region e = 1 and 



;: ~t~. \ x ~n) - o (m) 



y = nip (H4) 
x = spar (115) 



and in the second region e = — 1 and 



We can transform (113) into 



y = (n + 2)<p (H6) 
x = Vn~T2r. (117) 



x ^y + ^. +£ £ Xe v = Q (us) 

ax" 1 dx 



which is the equation considered in [6] for ee = 1 . 

Following H. Lemke we write down a solution to Eq. (118) in a compact form in 
three cases (concerning an integration constant introduced by H. Lemke). We adopt 
his solutions to our problem. 

1) C = k 2 > 0, k — arbitrary positive number: 

^ = - 21 "(f(©" + (l)")) +ln(4K2s2) (119) 

where a > is an arbitrary constant. 

2) C = —u) 2 < 0, oo is an arbitrary positive number: 

y(x) = -21n(2xsin(wlna; + 5)) +ln(4w 2 ) (120) 

where 5 is an arbitrary constant. 

3) C = 0: 

y(x) = -21n(|-ln(!-)) -21na (121) 
where a > is an arbitrary constant. 

Using these solutions we write down a spatial dependence of G e ff in the case of 
small and large fields. In the case of small fields one gets 

1) G . ff = GN (W) il±^i^ (122) 
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where 

r=^z\fnX^ p. (123) 

2) GtB = G ^ZZ n ^ 2( " +2)/ "( S in(^nO) 2( " +2)/ " (124) 

where 



2d 

ujlnS = -S (126) 



3) G e f f = G N a 2( - n+2 ^ n r 2( - n+2 ^ n (lnr) 2( - n+2 ^ n (127) 

and r is given by Eq. (123). 

In the case of large field one gets: 

1) G eS = G N (iK 2 a 2 ) { -±—f- (128) 

where 



r = —^(n + 2)X c0 p. (129) 

2) G eff = GwQ r 2 (sin(wlnr)) 2 (130) 

where 



r = -=A/(n + 2)A c0 p, In 5 = . (131) 

2d oj 

3) G e ff = G N a 2 r 2 (lnr) 2 (132) 

and r is given by Eq. (129). 

Let us notice that in that spatial dependence for large and small field we have 
k and a (5, uj) as integration constants. In this way integration constants induce a 
power law of this dependence and also a scale. For sufficiently big n (n > 14) there is 
no significant difference between both cases. It means the quintessence field behaves 
everywhere as for large field case (in these solutions of course). It is evident that the 
spatial dependence (in cylindrical symmetry case) of G e ff goes to some kind of the fifth 
force. However, we have to do not with a universal law of Nature but rather with some 
kind of initial conditions. Thus the real dependence of G e g on spatial coordinates can 
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be obtained after averaging G e s with respect to initial conditions. Let us describe it 
using Eq. (128). 

Let us write G e s in a form where k and a are explicitly visible: 

i+<:#y 
(t) : 

U(n + 2)W (134) 



G eff = G N (4k 2 a 2 ) • v / 2(K _ 1} y , (133) 



X 2 



The constants a and k are chosen randomly for they are dependent on initial conditions. 
Let p be a measure defined on (0, +oo) 2 , positive and normalized to 1, i.e. 



J J dfi(a,n) = 1. (135) 





This measure gives an account how frequently we have to do with some initial conditions 
(i.e. with integration constants a, k). In this way an experimental dependence of G e s 
on x should be obtained from 

+oo +oo K 2_ 2 / 1 + ( ^2^ 2 



E(G eS ) =AG N J J , > dp. (136) 

V a I 

It means it is an expectation value of G e ff with respect to the measure p. p need not 
be absolutely continuous with respect to the Lebesgue measure on M 2 . The formula 
(136) could give an account on some serious problems with comparisons of several 
measurements of G (a Newton constant). Maybe we have to do with different initial 
conditions for quintessence field. These initial conditions appear with different proba- 
bilities according to the measure dp(a,K) going to an expectation value E(G e s). The 
second central moment of dp(a, k) (if it exists) can express a deviation from the law 
described by E(G e s). 

In the simplest case we can suppose a continuous Gaussian distribution for k only: 

dpin) = -jL- e -(— o) 2 /(2<x 2 ) dK _ 



In this case we have 



-2 + °° (l + (^r) 2K Y 

G eS (x) = E(GMx)) = %^ / dKe-l"-"* > 2 /^V L ; 2 ;\/ . (136a) 

— oo V a I 
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One gets after some algebra 

G eS (x) = G n k exp ( a - \ y— J 



(2ln(^)a 2 + K 2 ) 2 /n 
4cj 4 exp ^ ^- In (-) 



2a 2 



+ (2a 2 + Kq) exp 



(21n(f)^ + ^) : 



(136b) 



2a 2 



4 (fJ + K ) 2 2 / X \ 2 2 A 

4cr exp- — In I — J + 2(a + Kq) exp 

Zi \ Ob / 



(a + k ) 5 



In the case of the solution (107) parametrized by two functions of one variable we 
can consider them as random variables parametrized by z and t. Moreover in this case 
we should consider a measure [i on an infinite-dimensional space of functions / and g, 
supposing that solution (107) is generalized to this space. Afterwards we can use as [i 
the Wiener or Gaussian measure in L 2 space. If we use for Gaussian distribution the 
normalized distribution N(0, 1), the formula (136b) simplifies to 



GW(x) = e n ^eG 



N 

^ xx 2(1 + Iu(,,,m) 

a 



ta (f) + (f) 



2(l+ln(a;/o)) 



-4 In' 



(f) 



+ 2 



(136c) 



Let us consider three our cases (119), (120) and (121) for small and large fields 
cases. 

The small field case is such that 



e v < 1, i.e. <p < 0. 



The large field case is if 



e v > 1, i.e. (p > 0. 
For (119) we have for the small case 

f(r) > 1, 

where 

/(r) =r 1 " K + r K+1 , r = (f) 
Let us consider two cases 



(137) 
(138) 



(139) 



< k < 1 and k > 1. 
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In the first case f(r) / in [0, +00) and we have simply r > r where r satisfies the 
equation 

r 2 K+1 + r - 1 = 0. (140) 

In the second case k > 1, 

lim f(r) = +00 and lim f(r) = +00. 

r— >0 r— >oo 



The function /(r) has a minimum at 



1 \ V« 
K — 1 



, \ k—1/k / -, \ k + 1/k 

K + l\ / « — 1 



Let us calculate f(ri). 



^ = H^ttJ ■ (142) 

It is easy to see that if k > 1, then k — 1/k > 0. This means that 

/(n) > 1 (143) 

and therefore /(r) > 1 for all r > 0. Thus simultaneously we get a solution for large 
field only if k < 1, i.e. 

r < r . (144) 

If k = 1, we always have /(r) > 1 (i.e. only a small field). 

Let us consider (120). In this case the small field condition reads 

x 

h(r) = rsin(ojlnr) > 1, where r = — . (145) 

5 

First of all we need h(r) > 0. Let us observe that \h(r)\ < r for every r > 0. Next we 
see that the roots of h(r) are the numbers 

r ,k = e k7T/uJ , fc = 0,±l,±2,... (146a) 

and 

h(r)>0 if r ,2fc < r < r ,2fc+i, fc = 0,±1,±2, . . . (146b) 
The maximum of the function h(r) in the interval (ro,2fc, ^o,2fc+i) is greater than 

Kr ,2 k+ i/2) = e 2k ^-e^\ (147) 
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but smaller than r^^k+i- Thus the maxima are smaller than 1 for negative k and greater 
than 1 if k = 0, 1, 2, ... . It means that in each interval (ro,2fc, r o,2k+i), k = 0, 1, 2, ... , 
there exist two numbers r 3t2 k and r 2j 2k such that r 3j2 k < ^2,2fc, 

/i(r 3 ,2fc) = /i(r 2 ,2fc) = 1 (148) 

and 

h(r) > 1 if r 3i2fc < r < r 2>2 fc . (149) 
The condition for large fields, < h(r) < 1, is satisfied if 

r ,2k < r < r ,2fc+i, fc = -l,-2, ... (150) 

or 

r ,2k < r < r 3j2 k, or r 2>2 fc < r < r 0>2 fc+i, fc = 0,1,2,... (151) 
In the third case, i.e. Eq. (121), one has for the small field case 

rlnr>l (152) 

where r = -=-. Let r 4 lnr 4 = 1, r 4 > 1. We have r > r 4 = 1.7632 .... In the large field 
case 

< r < r 4 = 1.7632... . (153) 
Thus we have in general large fields on large distances. 

In this way we have solutions for large and small distances. One can try to connect 
them to get a solution for all distances. However in this case it is necessary to be very 
careful, for our solutions depend on some integration constants which can be different 
for both asymptotic regions. 

Let us consider Eq. (63) in two special cases: 
I Qo = q(z) — static and depending only on z; 
II qo = qo(t) — non-static and spatially constant. 

Let us consider also these cases for small and large fields ip. In all of these cases 
we come to the following equation 

S + eie&v = (154) 

where e\ = 1 for case I and £\ = — 1 for case II. 
Eq. (154) can easily be reduced to the integral 



x _ Xn = £i [ d y n 551 

V2J Je 2 uj 2 - riev 



where r\ = e\EE, rf = 1, C = 2e2^ 2 , w > 0, e\ = 1 is an integration constant, e\ = 1 
and xq also is an integration constant. 
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A. 



B. 



After some calculation we get the following solutions: 

y/2 (x — Xq)lo 



y(x) = 2 
y(x) = 2 



In to — In 



sinh 



In to — In 



sinh 



\/2 (x — xq)lo 



, rj = -1, £ 2 = 1 
, f] = 1, £2 = 1 



where 



or 



^(*-*°)" > ln(l + ^2) 



>/2(- 



X — Xq)UJ 



< \n(V2- 1). 



(156) 
(157) 

(158) 
(159) 



2/0*0 = 2 



In w — In 



cos 



V2(x- 



x )u 



, 77 = -1, £2 = "I- 



(160) 



Let us apply these solutions to our problems. First of all let us consider a static 
configuration with z dependence only. In this case £\ = 1 and 7/ = ££. For the small 
field case one gets, e = 1, 77 = e, 



where 



P 



G 



N 



L0 2(n+2)/n 



|sinhp| 



(n+2)/n 



nA c0 . 



For the large field case we get, e = — 1, 77 = — £, 



G e ff = — ^ |cosp| 



where 



P= — A ° ; V(n + 2)A e0 . 



(161) 
(162) 

(163) 
(164) 



In this case e = 1 for 7/ = —1. 

In a nonstatic configuration E\ = — 1 and 7/ = — ee. For the small field case (e = 1), 
77 = -£, 



where 



a; 2(n+2)/n M 1 ?! 

e 3 u;(t - t ) 



(n+2)/n 



71 A, 



cO 



(165) 

(166) 
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and analogically for the large field case (e = — 1), rj = e, 

GW = — y \cosq\ , (167) 

q =^^4^j^. (168) 

In this case e = — 1 for r] = —1. 

Let us notice that in a static configuration for small field we have two possibilities 
for i] = —1 (no condition on p) and rj = 1 (conditions (158-159)). Thus without 
conditions we have e = — 1 and with conditions e = 1. In a non-static configuration 
for small field we have vice versa e = 1 without conditions and e = — 1 with conditions 
(158-159). 

Let us come back to the Eq. (89) and consider it in a travelling wave scheme. In 
this way we have 

q (z,t) =q(z-vt) (169) 

where v is a velocity of the travelling wave (a soliton), \v\ < 1. Let us consider this 

equation in both regimes (for small and large fields). In this way we come to the 
expression 

(1 - v 2 ) - eee x = (170) 



2 - eieee x = (171) 



where \ ls a shape function of a soliton. Changing an independent variable from £ to A 
one gets 

dX 

where £\ = —1, i.e. we get Eq. (154) with r] = —ee, 

A = - 7 =L=, £= Vl-« 2 A. (172) 

In this way we adopt our solutions A, B, C in both regimes: small and large field 
(changing % into A). For small field we get (e = 1, rj = —e) 

2 

qo(z, t) = =r- [In a; — In |sinhp|] (173) 

/3n 

where 



aVnAgo 

p =i7 _ (,-,*). (174) 



For large field we get (e = — 1, rj = e) 



2 

9oM) = =5; — — - [lnw-ln|cosp|] (175) 
p(n + 2) 
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where 



V 



u\J(n + 2)\ 



(z - vt). 



(176) 



For (173) we have rj = —e and because of this e = 1 without any conditions and if 
e = — 1 we have conditions (158-159). In the case of the formula (174) r\ = — 1 and 
e = -l. 

We can write down formulas for G e g in the soliton case 



G e ff = 



G 



N 



L0 2(n+2)/n 



|sinhp| 



(n+2)/n 



(177) 



and p is given by the formula (174) (with or without conditions (158-159)). This is of 
course a small field case. 
In the large field case 

G eff = ^ \cosp\ , (178) 

and p is given by the formula (176). In this case e = — 1. (Let us notice that this is a 
case of SO(3) group in our theory.) 

Let us notice that conditions (158-159) can be considered as conditions for small 
field in z or t domains. Let us notice that in our solutions concerning a behaviour 
of an effective gravitational constant we get completely arbitrary length or time scale 
(given by integration constants). In this way a spatial or time dependence of G e ff can 
be (except the solution (80) and simultaneously the approximate solution in the case 
of spherical symmetry) such that G e ff can be really constant on distances (or times) 
accessible in experiments. Only a statistical approach mentioned here can give a light 
on this dependence to compare it with an experiment. 

In order to connect our results to the ordinary gravitational physics we consider 
again Eq. (67) in small field regime for initial conditions (p(0) = and (0) = 0. 
The first condition means that we want to have G e ff(0) = Gn and the second that the 
quintessence field does not grow quickly. The problem cannot be solved analytically. 
Moreover R. Emden in the first point of Ref. [5] did it for us. We quote here his results 
adopted to our notation (e = +1 and e = —1). 



X 


-y 


ey 




0.00 


0.00000 


1.00000 


1.00000 


0.25 


0.01037 


0.98969 


0.98823 


0.50 


0.04113 


0.95971 


0.95409 


0.75 


0.09113 


0.91290 


0.90109 


1.00 


0.15903 


0.85296 


0.83380 


1.25 


0.24225 


0.78486 


0.75816 


1.50 


0.33847 


0.71285 


0.67920 


1.75 


0.44488 


0.64090 


0.60143 
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X 


—y 


e y 


ell / i V 


2.00 


0.55967 


0.57140 


0.52749 


2.50 


0.80584 


0.44671 


0.39813 


3.00 


1.06226 


0.34537 


0.29670 


3.50 


1.31937 


0.26730 


0.22138 


4.00 


1.57071 


0.20790 


0.16611 


4.50 


1.81246 


0.16325 


0.12601 


5.00 


2.04264 


0.12968 


0.09686 


6.00 


2.46598 


0.08493 


0.05971 


7.00 


2.84160 


0.05833 


0.03887 


8.00 


3.17489 


0.04180 


0.02656 


9.00 


3.47128 


0.03108 


0.01893 


10.00 


3.73646 


0.02384 


0.01398 


100 


8.59506 


0.000175 


5.0854 • 10" 5 


1000 


13.09847 


0.000002 


3.0683 • 10- Y 



where 



x = \ \/ nX cor 



lOMpc 



g eff = ( e s/)("+ 2 )/™ = ( e v) 8 / 7 
G 



N 



(179) 
(180) 



We take n = 14. 

It is easy to see that for large n 



G e ff 

G~ 



N 



(181) 



It is easy to see that on a distance of 1 Mpc G e s does not differ from Gn- Even 
on a distance of 10 Mpc it is about 10% smaller. Thus in the Solar System Newtonian 
gravitational physics does not change. Even on the level of a galaxy this change is 
minimal and cannot be observed. Moreover, there is an important conclusion: on 
distances about 200 Mpc the strength of gravitational interactions is about 10 -5 times 
this on short distances measured in the Solar System (and for 10 3 Mpc of 10~ 7 ). It is 
hard to tell how it influences a mass of a cluster of galaxies if we realize that from any 
observational data only a product GM has been obtained (not M). 

From the other side on distances of 100 Mpc the strength of gravitational interac- 
tions is very weak (not only because of the distance). Thus if we consider clusters of 
galaxies as substrat particles in cosmology then they do not interact. 

Let us consider Eq. (63) in Cartesian coordinates for two-dimensional static case 
(i.e. £ = 0, §- t = 0). One gets 



d 2 q d 2 q 



dx 2 



dy 2 



ea exp(n/3(7o) (exp(2/3g ) — l) = 



(182) 
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(where a, [3 are given by formulas (64) and (65)). 
As usual, we come to the formula 



(S+i)-^- 1 )- 



where 



'X, 



i = l,2. 



cO 



We consider Eq. (183) for small and large fields and we get 

'd 2 x , d 2 x\ ~ x _ 
— eee x = 



dz\ dz\ 
where as usual for a small field e = 1 and 



and for a large field e = — 1 and 



X = nip, 



X = (n + 2)ip, 
Zi = y/n + 2xi 



Thus we come to the equation known as Liouville equation 

Ax = e x 



(183) 



(183a) 



(184) 



(185) 
(186) 



(187) 
(188) 



(184a) 



if ee = 1. 

This equation can be explicitly solved. First of all we change independent variables 

into 



Z = —(z 1 + iz 2 ) 
x(Z) = -ln(I(l-| 5 | 2 )) + iln 



dg 



dZ 



(189) 
(190) 



where g is an arbitrary analytic function on a complex plane Z . 
In this way we get for the small field case 



G e ,ff — G 



N 



2X (n+2)/n 



dZ K ' 



-(n+2)/(2ra) 



(191) 
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where 

In the large field case 

G e ff = Gn 

where 




— ■(x + iy). 



A 



cO 



2(n + 2) 



(x + iy). 



-1/2 



(192) 



(193) 



(194) 



Eqs (191) and (193) can have very interesting behaviour for could have some sin- 
gularities. The physical interpretation of these singularities can be very interesting. 

Similarly as for Eqs (104) and (107) we can consider an expectation value of G e s 
with respect to some kind of normalized measure for a space of analytic functions on 
the complex plane (respectively chosen). 

It seems that an assumption all the energy of a quintessence is stored as quint- 
essence particles is unrealistic. Let us suppose that only a fraction of this energy is 
stored as particles. Let this fraction be 77, 



< 77 < 1. 



(195) 



r\ can be a function of time and even of a space-point (locally). 

Let us suppose that a gas of quintessence particles is a perfect gas governed by 
the Clapeyron equation. Thus we have 



Pi 
Pi 



K B T 

m 



T 



T n = 



m 
K B 



0.11°K 



(196) 



(if we take mo — 10 5 eV). T is the temperature of a gas. 
One gets 



Pi 
P 

Eventually one gets 



VPQ = VP 
PQ+Pi = 



T 

Pq + ri — p 
Jo 



-p(l -ri)+p — r). 



P 
P 
w 



= P [V [1 + 

wp 



T 



- 1 



(197) 
(198) 

(199) 
(200) 
(201) 
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Now we can calculate an isothermic speed of sound in a quintessence. 



cl 



p 
P 



w 



V 1 + 



To 



- 1. 



(202) 



For 



we get 



< C\ < 1 



T 



T + T 



< rj < min 



2Tn 



T +T' 



(203) 



(204) 



Let us remind to the reader that an isothermic sound is appropriate for low frequency 
of acoustic waves (we have to do with this sound in astrophysics) . In general we have to 
do with so called adiabatic sound. In order to calculate a speed of an adiabatic sound 
in a quintessence we should find an analogue for a Poisson adiabate for our equation 
of state. One gets supposing that an internal energy of a quintessence is an energy of 
one-atomic gas of quintessence particles 



dU+pdV = 



(205) 



or 



and finally 



— = const. (207) 



3 

and a speed of an adiabatic sound simply reads 



K = 5L±1 (208) 



c '=y?=f G( i+ £H- <2o9) 

It is interesting to ask what kind of a polythrope k represents. Let us remind to the 
reader that in general 

k = (210) 

c v c 

where c is a specific heat of a polythrope in mind. 
One gets in our case 

c = 3( r 1} r (2ii) 

2r/ v ; 
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where R is a universal gas constant. 

Thus we found both speeds of sound in a quintessence for low frequency and for 
high frequency of sound. The measurement of both speeds can help us to find 77 and 
T. One gets 

'-!((§)"-) 

and 

T = 3(cq- ci) ( cl(2Cl + 5) " 3C|) • (213) 

Let us notice that a gas of quintessence particles is quite cold. Moreover these particles 
are highly relativistic. For the temperature T~T = 0.11°K one sees that a speed of 
a quintessence particle is about 0.91 of the speed of light. Moreover we can consider 
lower temperatures. It is interesting to notice that the speed of sound is of the same 
order. 

Thus if we want to have both speeds smaller than a speed of light, 

0<Cf<l, i = 1,2, (214) 

we get 



where 



The condition 



1 + V13 + 12t 

v > — JT^ = f ® ( } 

t = Y Q - ( 216 ) 

t = y > 13 (217) 

guarantees that f(t) < 1. Moreover for t = 14 one gets 

77 > 0.96357. (218) 

This seems quite interesting, however maybe too much. Let us calculate a mean scat- 
tering length of quintessence particles for such a big rj = 0.96357 

'scattering = = 1-0378 ■ I. (219) 

rjan 

We have still to do with a very dense gas of quintessence particles. 

It seems that it would be reasonable to repeat these calculations using a different 
equation of state for quintessence particles gas. In particular we can consider a gas 
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of quintessence particles as a massless boson gas (with spin zero). In this case the 
equation of state looks: 

P = PQ + Pi = "(1 - V)P + T 4 (220) 
and an adiabate equation 

dU = c v dT + P dV = (221) 

where 

c v = 16crT 3 , (222) 



a is the Stefan-Boltzmann constant. 
One easily integrates 



■4 = const. (223) 



where 

(224) 

We have as before two kinds of sound: an isothermic sound 



r '=v! = v''- i+ T T4 (225) 



and an adiabatic sound 



B, = ,/»(i-i + ^r«)=^- P») 



2 

From < Cj < 1, i = 1, 2, one gets 



1 

— — j < 7/ < mm 



1 + (T/T ) 

and 



i, 2 



7^ \4 



i + (t/t ) 



(227) 



"> " 3t4 " 1+ 2(f + ^) 18t4 --7w ( 22S ) 

with the condition 

T 1 

t = > —= ~ 0.707. (229) 
T y2 

This condition guarantees that 

lit) < 1 (230) 
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where Tq = ^ or 



Kb is the Boltzmann constant. 

From the measurements of Ci and C 2 we can obtain as before r] and T. One gets: 



(232) 




and 



Recently many papers have appeared concerning the speed of sound in a quintessence 
(see Ref. [7]). In some of them the authors propose to measure this speed. 

There are interesting propositions to include primordial gravitational waves to 
fluctuations of a quintessence and vice versa. We examine these problems below. 

We consider a mass of a quintessence particle, a speed of sound in a quintessence 
and several solutions to quintessence equations coming to the interesting behaviour 
of an effective gravitational constant. Some of issues of our theory are considered in 
self- interacting Brans-Dike theory [8]. 

The fraction r] can be connected with the part of an energy density of quint- 
essence field Q in such a way that it is a fraction of fluctuation energy density around 
an equilibrium Q . In this way the qo field which can evolve due to an evolution of Q 
and due to fluctuations of primordial gravitational waves will be a source of a gas of 
thermalized particles. It seems that an approach with a boson equation of state is 
more appropriate to consider. 

It is interesting to consider a fraction of a dark energy stored as boson particles 
as a dark matter in the Universe. 

Let us consider a primordial spectrum of gravitational waves in our approach. 
This spectrum is flat for our Hubble constant is really constant 

Pgrav = ( 2^) (235) 

(see Ref. [9]). Thus 

n gr = 1. (236) 

Now we start to examine quintessence fluctuations caused by fluctuations of a 
metric. In order to do it we follow Ref. [9] to perturb a spatial part of a metric 



R 2 (t) -dr 2 + (Sij + 2Ej^ dx l dx j , (237) 
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r is a conformal time. 

If we expand Einstein equations for (237) up to linear terms, we get 

d 2 E T t dE T t 

, l jfc + 2RH K ^± + k 2 E T t = 0, if = 0,1, (238) 
dr 2 dr 

where Ej- is a spatial perturbation of the metric, R{t) is a scale factor depending on 
a conformal time, E T . t is a k Fourier component Ej- 

V ( 2vr )3/2 J %3 k ^ 239 ^ 

k = (&!, k 2 , k 3 ), \k\ 2 = k 2 , f=(x,y,z). 

The important quantity is an amplitude of gravitational wave corresponding to Ej-, 
i.e. 

h..r = RE T r . (240) 

V k ijk y ' 

Hk is a Hubble constant. We take K = 1. The relation between a conformal time r 
and ordinary time t is given by 

i2(r) = --!-, -oo < r < 0, (241) 
i?(i) = i? ie fflt . (242) 

For /i . . r one gets 

2 d/t.r , 

_!jfc +fc 2^_ 243 

This equation can be easily solved. For further investigations we need only 

h s = Kk ( 244 ) 

h % = A % (r\k\ cos(r|fc|) - sin(r|fc|)) + B % (r\k\ sin(r|fc|) + cos(r|fc|)) (245) 
and especially 

^ = ^rsin(r|/c|) +^rcos(r|A:|) (246) 

where and i?£ are constants. 

We need also a quintessence field time dependence (e.g. in a slow roll approxima- 
tion). We use our solution from the last point of Ref. [2] (see Eq. (14.416)), making 
some simplifications and changing t into r: 
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Moreover what we really need is a derivative of 



l-C 2 (-r) 2 « 1 J 



where 



k = —j- — (248a) 
Hi 

C = e *^ t0 {RiHif^'" 1 = e ^v^*o (R x H{f . (248b) 
Now we proceed to a quintessence fluctuation equation (see Ref. [7]): 

5Q % + m^Qj: + (c 2 k 2 + R 2 U"(Q)) 5Q % = ^=h' k ■ V . (249) 



One gets 



5Q % + SHtfQz + {c 2 s k 2 + R 2 U"(Q)) 5Q U 

Ck f (-t)*- 1 sin(r|A;|) (-t) s " 1 cos(t 

= -=r- Ar — : h Br 



V k l-C 2 (-r) 2k k 1-C 2 {-t) 2R J (250) 

CR 8in(r|fc| + jg) 2R 
/3 fc l-C 2 (-r) 2S 1 J ' 

/3 is a normalization constant in a definition of Q = =, i = 8 ^|^ (see (28)). 

Let us consider 5Qr as a Fourier component of a go field subject to the following 
initial conditions: 



^0fc(°) 

dr 



= Q O k(°)- ( 251 ) 
In this case one gets 



c s is a speed of sound in a quintessence, Cg and 5^ are constants. Let us notice that 
r = corresponds to t — > oo. Thus in some sense we have asymptotic conditions for 
gravitational waves. 

For further convenience it is good to change r into — r. In this way we get 

rG(0,+oo), r = -i-e- Hl *. (254) 
itl-Hl 
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The fluctuations (the field qo) of a quintessence Q are driven by gravitational 
waves and a quintessence field Q. It is easy to see that we are interested in solutions 
for small r (large t). In this case one gets 

^f- 3H ^ + hJT^ = ~^P C ^\k\ - 5,). (255) 
Let us calculate a constant k. One gets 

B^hi ( n \ n/ Vl7lV /4 v^ 



K = 



4(n + 2)| 7 | 1 / 2 



Hi \n + 2 / \ (3 J rn p i 

n \ n/4 /m-\ (n+2)/2 IPIV2 / I Si \ n/2 (256) 
= 8v^(n + 2)V / M(-^j) ( ^ ) ^TY " ■ 

If we use our simplified model with n = 14, a 2 . = a em , = rasw, we get 

k^IO 92 . (257) 

Thus this is really a large number. 
Moreover for small r we have 



In this way we arrive to an equation 



and finally 



g - (T ) = - _ B ^ CK r *+ 2 - _ A * CK t* +1 (260) 
ofcV 7 /3(k + 2)(k + 1) P(k + 1)k 



or (taking under consideration that k is very large 

C 



(r) = -^(B s fcr s + A £ r s - 1 ) (261) 

+ fe^t^ff!)-" 1 ^/^'- 1 ^- 1 ^ • (262) 
g (i,r) = 1 q^{ty %F dk. (263) 
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One gets 



<lo(t,r) = -^ g(r) e~ B ^ 



(3k \ \Ri 



-"1 \ _-2BV2^lt c H x t 



(264) 



+ m ( ^ ) e-™«™ ■ e 



where 

ikf j3 



f(r) = J A^ r d A k (265) 

g(f) = J {klB^ d 3 k. (266) 

g(f) and f(r) characterize a spatial dependence of gravitational waves background. 
For large t we have to do with g(f), i.e. 

go ( t> ^ = -£L g (f)( RlHl )-(@y&;/H 1 )-i) . e -Bv^t e H lt (267) 

(3k 

In this way one gets from an energy-momentum tensor for q$ 

Tpv = d^qo • d„q - ^rj^ (d a q ■ d a q + mlql) (268) 

P qo = T44 = - - (ql ~ | Vgol 2 - mggg) 

= \ql ( (Hi - 5^) 2 + | V In g(r) | 2 + mg) (26g) 



p sv 2nB 2 Rj ((Fx - BV2^) 2 + | V ln 5 (r)| 2 + m§) 



2 p 2BV2^t 



(3 2 e 



where 



V (270) 



is an energy density of primordial gravitational waves. 
For k = B fl^™ is a large number, 

BV^l > Hi (271) 

and p qo is going to zero if t — ► oo. This density will be frozen (no dependence on time) 
at an end of the second de Sitter phase t = t]} nd , to = initial- ^ n * n ^ s wa y one S e ^ s 

2nS 2 J R?((F 1 -5 v ^) 2 + |Vln 5 (r)| 2 + m2) 

p gp V_ L e 2Ar i (272) 

p ST ^exp( 2 ^ jVl ) 
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where Aq is an amount of inflation during the second de Sitter phase. If g(r) = const., 
p qo is isotropic and homogeneous. 

The function g{r) can be written in a form 

g(r)=go + Sg(f) (273) 

where Sg(r) is a small deviation and go = const. One gets 

V (In g(r)) ~ — V5g(f). (274) 

In this way one writes 

Pqo = Pq + Spqo (275) 

where 

2nB 2 Rlp sr ((#1 - BV2^) 2 + mg) ^ 
Pq °~ ^ex P [^^) 6 

= | V <?(f)| 2 2ni fg^ e 2^ 
^ex P (Mypv 1 ) 

According to our ideas the energy p qo should be stored as a gas of quintessence 
particles. The measurement of a low frequency speed of sound and a high frequency 
speed of sound can give us r\ and its spatial variation. In this way we get also a square 
length of a gradient of Sg(r). If we suppose that p gr is stored as a gas of gravitons 
(massless particles with two polarization states) then one can write 

Pgr = Pgr (^-) (277) 



(276) 



(see Eq. (231)). 

Probably we should also consider a gas of skewons. However, this is a different 
story. Taking into account (235) and (236) we should expect 5p qa as extremely small 
and we can neglect it in our theory. Moreover, /> gr is very important. Let us consider 
Eq. (253) for these values of r for which the variable t is close to t* nd (the end of the 
first de Sitter phase or beginning of the second de Sitter phase) . In this way one writes 

e -HitL d + £ = To + £ 0<e«l- (278) 



H\Ri 

Taking into account that £ is very small we arrive to the following equation 



3 Hl + ^ + ^ = -||sin(^| +? s ) (279) 
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where 

h = ^\ k \~h- ( 28 °) 

We have also initial conditions 



d %k 



(°) = = °- ( 281 ) 



Eq. (279) can be solved by a Laplace transform method in the case of 5^ = 0. 
Let 

m, 

U2 2 
n l T 



-,2 



W + ^ = b 2 . (282) 



One gets 



2/3 V * 



a 2 9 0j f(a) - SH l8 q on (8) + &\ g (s) = arctg ( U. ) (283) 



and 



a As) ~ arCtg ^) (284) 

C c r /" arctg(^)e sS , 
o 

However, we cannot get any compact form of this solution. Moreover for we are looking 
for a solution around zero, for 5^ = we get 

-°& Mma -°S!M. (286) 
In this case the Laplace transform method is successful: 

OO 

-3^/ , (287) 

o 

where s\ and S2 are roots of the polynomial 

s 2 -3H lS + b 2 = 0. (288) 
Moreover £ is small. Thus we write 

e^^l + s^, i = l,2, (289) 
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and get 



if 



A = ml - Ab 2 = ml - i(c 2 s k 2 + -gO > 0. (291) 
If A < 0, one gets 

%d0 = ~^=f (l + eMl^O (^shiK) _ cosK))) (292) 

where 

4^ 2 = -A = A(c 2 s k 2 + - 9Hl > 0. (293) 

For small £ one gets 

go£(0 = - ^ - e, (294) 
the dependence on i is as follows: 

e = — ^(e-^-e-^ nd ). (295) 
Moreover t is close to i* nd and we get 

^— (*-4d) (296) 

for i very close to ig nd . 

In this way we get fluctuations of a quintessence in the moment of the first order 
phase transition of the first de Sitter phase to the second de Sitter phase. 

Finally one gets 

_ _ C # |e-^ 



X 1 + 



m x 



(297) 



^9# 2 - 4( C 2|fc| 2 + i22mge 2f?1 *end) 



for 



|fc| < —\Jm 2 - m\e 2H <^ml = k (298) 
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with a condition 



3ffi > 2i? 1 e Hl *- d m 



(299) 



and 



One finds 



e ~^ 1 *end 



( 



1 



3#! 



v J9^f - 4 (c 2 ^ 2 + i22m§e 2Hl *end) ) 



d 3 ke ikf! 



\k\>k 



(c 2 s k 2 + Rjmy^l^ 



(300) 



(301) 



Equation (301) gives us space-time fluctuations of a quintessence field after a 

phase transition from the first to the second de Sitter phase. Let us notice that (297) 

is singular for k — ► /co- 
in an approach developed here we consider Eq. (250) with zero initial conditions 

in various approximations. Moreover there is a different approach. In this approach 

we write a solution to Eq. (250) 



9 fc(*)=9ofc(*) + 9oiE(*) ( 302 ) 

where %%{t) is a solution of a homogeneous equation (a general integral) and %%{t) is a 
special solution to inhomogeneous equation. This different approach consists in taking 
%k (*) = and considering only q ^(t). In this case we get nonzero initial conditions for 
quintessence fluctuations. However, this approach is not unambigous. Our approach 
is unambigous. The zero initial conditions are in some sense distinguished. 

Let us calculate an energy density of this scalar field. In order to do this let us 
write qo(f,t) in a form 

?o(f, t) = A(r) + B(r)(tl nd - t) (303) 



39 



where 



A(r) 



B(r) 



WRJh J 



3 I. „ifcr 



cPke 



|fc|<fc 



1 + 



cm 

3iJi 



V 



9if? - 4 ( C 2fc 2 + R 2 m 2 e 2H <^ 



3 e -^l<end 
4/3i? 



end f 

r 7 



.fcr 



|fc|>fc 



(c 2 k 2 + R 2 m 2 e 2H ^^ 



One gets from the Eq. (268) 



(304) 



(305) 



where 



P qo = Tm = g (C(r) + £>(fO(Ci - <) + E(mL d - t) 2 ) 



C(r) = B 2 (f)+m 2 A 2 (f) + | V(A(f))| 2 

D{r) = 2 ( V (A(r)) • V (5(f)) + r^A(r)S(r)) 

£(r) = (|V(S(r))| 2 + m^ 2 (r)) . 



(306) 

(307) 
(308) 

(309) 



All these formulae are satisfied for t close to t* nd . In general p qo is anisotropic and 
inhomogeneous. For t* nd — t is very small we get 



(310) 



It is interesting to calculate a time t(r) for which p go is zero (t(r) depends on a 
space point). One gets (neglecting gradients of A(r) and B(r)) 



<<' r) = t '- + 4( 7( 



r) + =- 



/(0 



(311) 



where 



3Hi 
2R~[ 



!\k\>k d3ke 



ikr 



3l 



( C 2fc2 + R2 m 2 e 2 ^l* end ) 



I\k\<k d3k ' 



.ikr 



C,\k\ 



( C 2fc 2 +R2 m 2 e 2H l* end) 



1 + 



ZH 1 



9H 2 -4(c 2 fc 2 + it 2 m 2 e 2 - fflt end) 



(312) 
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The conclusion from these calculations is that after a very short time an energy 
density of quintessence fluctuations is going to zero (almost zero) and probably is frozen 
on a very low level. A time of this process depends on a space point. Moreover, after 
sufficiently long time (which is really very short) we get a homogeneity and isotropy, as 
far as our rough calculations advise us. Let us come back to Eq. (253). Our solutions 
obtained here are for large t and t close to t J end . One can try to match them. It seems 
reasonable to suppose that p° (see Eq. (276)) is equal to this very low level energy 
density, mentioned above. 

The quintessence field qo(f,t) considered above can be a source of fluctuations of 
an effective gravitational constant G e ff. 

G eS = G N exp(-(n + 2)pq (r,t)) (313) 

(see Eq. (28)). If we use Eq. (267), we get 



n n ( y/2(n + 2)^/3 ' 

& e ff = &n exp _ = 

((#! - B^Mf + I V lng(r)\ + mg) . 



(314) 



and we connect contemporary fluctuations of an effective gravitational constant to an 
energy density of quintessence fluctuations. Taking a quintessence field from Eq. (303) 
we get 

G eff = G N exp (-(n + 2)0 {A(r) + B(r)(t l end -t))). (315) 

In this way an effective gravitational constant depends on a speed of a sound in a 
quintessence. 

We can consider some applications of a quintessence in astrophysics. First of all 
we can develop a formalism concerning a macroscopic flow of a quintessence gas (i.e., a 
relativistic hydrodynamics of a quintessence) . For a quintessence field has an influence 
on an effective gravitational constant, such a flow can disturb gravitational interactions 
between galaxies or even stars. In this approach we should add some assumptions, e.g. 
that a quintessence energy density is an energy density of quintessence field in statistical 
field theory. However, this is beyond a scope of this work. 

Secondly, it is interesting to consider a Bose-Einstein condensation of quintessence 
particles. They are massive scalar particles. Under some assumptions (beyond this the- 
ory) we can find a wave function of a condensat. This wave function can be considered a 
field of a quintessence and afterwards enters the formula for G e ff. In this way Einstein- 
Bose condensation can influence effective gravitational interactions between galaxies 
and stars. This is also beyond a scope of the work. 

As a typical astrophysical application we can consider a problem of N bodies 
with the quintessence field. In this problem N points interact gravitationally with an 
effective gravitational constant depending on a distance. The simplest case is of course 
a two bodies problem. However, this problem cannot be solved analytically. It can 
be considered numerically using some codes for the N bodies problem. This is also 
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beyond a scope of our work. Moreover, some numerical solutions can be applied for 
galaxies movement in a cluster of galaxies. 

In all of those problems we should consider quintessence fluctuations developed 
here as a source of quintessence gas and Bose-Einstein condensat. In the third problem 
we should also consider an interaction of a quintessence field with mass points. 

It is interesting to consider a continuous distribution of a matter (a dust) under 
gravitational interaction and a quintessence. In this way we consider hydrodynamic 
equation coupled to a quintessence field in a Newtonian physics limit (i.e. selfgravitating 
system with G e s depending on a quintessence). 



Conclusions 

In the paper we consider a mass of a quintessence particle, a speed of sound in a 
quintessence and several solutions to quintessence equations coming to the interesting 
behaviour of an effective gravitational constant. Some of issues of our theory are 
considered in self-interacting Brans-Dike theory [8]. 

The fraction rj can be connected with the part of an energy density of quint- 
essence field Q in such a way that it is a fraction of fluctuation energy density around 
an equilibrium Qq. In this way the qo field which can evolve due to an evolution of Q 
and due to fluctuations of primordial gravitational waves will be a source of a gas of 
thermalized particles. It seems that an approach with a boson equation of state is 
more appropriate to consider. 
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